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We discuss a pseudospin representation of the two-dimensional t-J model. We introduce pseu- 
dospins associated with empty sites, deriving a new representation of the t-J model that consists of 
local spins and spinless fermions. We show, within a mean-field approximation, that our represen- 
tation of t-J model corresponds to the isotropic antiferromagnetic Heisenberg model in an effective 
magnetic field. The strength and the direction of the effective field are determined by the hole 
doping 5 and the orientation of pseudospins associated with empty sites, respectively. We find that 
the staggered magnetization in the standard representation corresponds to the component of mag- 
netization perpendicular to the effective field in our pseudospin representation. Using a many-body 
Green's function method, we show that the staggered magnetization decreases with increasing hole 
doping S and disappears at S ~ 0.06 — 0.15 for t/J = 2 — 5. Our results are in good agreement with 
experiments and numerical calculations in contradistinction to usual mean-field methods. 

PACS numbers: 71.10.Pm, 71.27.+a, 74.72.-h, 75.10.Jm, 75.50.Ee. 



I. INTRODUCTION 

It is well known that the essential physical properties 
of the copper oxide materials are described by the two- 
dimensional t-J modeli^. The action of this model is 
restricted to the single occupancy sector of Hilbert space, 
X) ( j c L c io' — 1- At half-filling, the t-J model reduces to 
an antiferromagnetic Heisenberg model. In the doped 
case, on the other hand, due to the single occupancy con- 
straint, electrons can move only onto empty sites. Thus 
the constrained electron operators no longer obey the 
fcrmionic anti-commutation relations. 

The slave particle mean-field methods 3 -^^^ were 
the first approaches to overcome this problem. In 
these methods, the constrained electron operators are ex- 
pressed in terms of auxiliary fermions and bosons. How- 
ever, as the local single occupancy constraints are re- 
placed by an average global constraint, these mean-field 
methods lead to considerable errors: for example, the 
antiferromagnetic long-range order is overestimated^^. 

Another approac h 10 ! 11 ' 12 ! 13 ] 14 is to express the con- 
strained electron operators in terms of spinless fermions 
and spin-1/2 operators by introducing pseudospins state 
associated with empty sites. Using this representa- 
tion, the t-J Hamiltonian can be described without any 
constraint in contrast with the slave particle methods. 
The spin- "up" and spin- "down" symmetry (time-reversal 
symmetry) of this representation is discussed by Wang 
and Ricei 2 - and Loosi 3 -. In the present paper, we develop 
this idea and calculate the critical doping S c , where the 
antiferromagnetic long-range order disappears. 

Our representation of the t-J model corresponds to 
the isotropic antiferromagnetic Heisenberg model in an 
effective magnetic field within a mean-field approxima- 
tion. In order to calculate the staggered magnetization 
of the model, we employ the many-body Green's func- 
tion method developed by Frobrich and KuntzA^. Un- 
like the slave particle methods, we show that the criti- 
cal doping S c for the disappearance of antiferromagnetic 



long-range order is in good agreement with the numerical 
calculationsi^Li^^. 

This paper is organized as follows. In Sec. |TT1 we in- 
troduce a transformation mapping the original Hilbert 
space of the constrained electrons to the tensor-product 
space of the spinless fermion and spin states. We also dis- 
cuss the difference between our method and the preceding 
method a 12 ' 13 . In Sec. IIII1 using the many-body Green's 
function method, we calculate the staggered magnetiza- 
tion of the model obtained in Sec. [TT] within a mean-field 
approximation. A summary is presented in Sec. IIV1 



II. PSEUDOSPIN REPRESENTATION 

We consider the t-J Hamiltonian 

H = -tj2 (4o c i<? + h - c -) -t^jL°^ + J J2 St ■ b j'( 1 ) 



<IJ>CT 



<hJ> 



where Ci a is the electron annihilation operator with 
spin a at site i, fj, is the chemical potential, and = 

\ E a ^ c L r a^ c i/9 1S ^ ne s P m operator with the Pauli ma- 
trices r. The sum is taken over all nearest neighbor 
bonds. The electron in this model is subjected to the 
constraint of no double occupancy, i.e., 



< 1. 



(2) 



Consequently, the basis of possible states in the t-J model 
consist of the states \ia) and iO), which correspond to 
a site singly occupied by an electron with spin a and 
to an empty site, respectively. The constrained electron 
operators Ci„ act on this basis as 



,\ia) = \i0), 4|<Q) = \ia). 



(3) 
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A. Pseudospin associated with empty site 

One of the most popular techniques to handle the 
single occupancy constraint ^ is the slave particle 
method 3 -^^^. In the slave boson representation, the 
constrained electron operators are given by the mapping 
c ia —> b\fi a , where bi is the bosonic operator annihilating 
the empty state and fi a is the fermionic operator anni- 
hilating the single occupied state. In this case, the non- 
holonomic constraint {5} is replaced by the holonomic 

constraint b\bi + J2a ftrfia = 1- F° r P rac tical purposes, 
however, these local constraints for each site are almost 
inevitably replaced by an average global constraint, re- 
sulting in errors: one of the most serious drawbacks is 
that the antifcrromagnetic long-range order is overesti- 
mated^^. 

In this paper, we develop another approach, rewrit- 
ing the t-J model in terms of spinless fermionic opera- 
tors (charge degree of freedom) and local spin- A opera- 
tor o 1 ^ 11 ' 12 ! 13 ' 14 . In this approach, the Hilbert space is 
mapped onto a tensor product space, \-)h ® \-)s, as 

\ia) -» \iO) h \ia) s , \i0) - |»l) h |iS)s, (4) 

where 

|iS)s = Ci|iT>s + C 2 |ii) S) (5) 

with |Ci| 2 + |C 2 | 2 = 1, is the pseudospin state associated 
with an empty site. We introduce the fermionic "holon" 
operator hi and the "spin" operator Si acting as 



hl\iO) h = \il) h , hi\il) h = \iO} h , 
St\H)s = \iVs, Sr\iVs = \il)s. 



(6) 



Owing to Eq. (j4]) , the constrained electron operators are 
written as 



Ci ^^ = h\{C x StSr + C 2 Srl 
c il ^ck L = ht(C 1 S+ + C 2 SrS+). 



(7) 



Clearly, the new operator Ci a has the same action as 
the operator Ci„ on the basis vectors ((4]). The four- 
dimensional space \-)h O \-)s consists of the orthonor- 
mal vectors |j0)^|icr}5, \il)h\iS)s, and \il)h\iS)s with 
\iS)s = Ci\i i)s — C||it)s- However, it is easy to prove 
that the operator Ci a acts on the \il)h\iS)s state as 



c ia \il) h \iS) s = 0, clJiljhliS) 



0. 



(8) 



Consequently, the unphysical basis vector |il)^|iS')s can 
be ignored when the action of is considered 1 ^. 

In this pseudospin representation, the "true" spin op- 
erator and the constrained electron number operator are 
expressed as 



E4 



(l-^)S,., (9) 

£4gi, = l-hK (10) 
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FIG. 1: Bloch sphere describing the orientation of pseu- 
dospins. 



Consequently, the single occupancy constraint 

XXW < 1 ^ 1 - ^ < 1 en) 

is automatically satisfied. Thus no approximations for 
the constraint are necessary unlike the slave particle 
methods245&LS£. 

B. Pseudospin representation of t-J model 

In our pseudospin representation, the t-J model is 
found to be expressed as 

H = -tY,[hih){Cl{StSrSfSj + S~S+) 



<i,3> 



~\-C\G<i (S^~Sj~Sj + Si S^~ Sj~*) 
+CtC 2 (StSrSr + SrSrSf) 

+C 2 2 (SrS+Sr S + + S+S j )}+h.c. 
-nY^h\hi • ./^ //,//, S, • Sjhjh] 

i <i,j> 



(12) 



which can be rewritten in the following more convenient 
form with Ci = cos § and C* 2 = e lv sin § [Fig. Q] 



<i,3> 



n- | -Si + -Sj + iSi X Sj 



+ h.c. 



<hj> 



where n = (sin 9 cos <p, sin 8 sin ip, cos 9) is a unit vector 
in the orientation of pseudospins. Loos^ put a certain 
restriction on the orientation of the pseudospins so that 
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the time-reversal symmetry of the t-J model is preserved, 
assuming that the "spin" operators Si are transformed in 
the time-reversal operation similarly to the "true" spin 
operators s,;. In contrast, we put no restriction on the 
parameters 9 and ip. In Sec. IIII B[ we show that our 
results are independent of the parameters 9 and ip. 

The representation of the t-J model introduced in 
Ref. [ll] can be obtained by setting 9 — in Eq. (fT5)) . 
The similar expression is also introduced for the Hubbard 
model in Ref. [lj. For simplicity, in this paper, we also 
choose the parameters as 9 = 0, i.e., the pseudospins are 
aligned along the z axis. Consequently, the pseudospin 
representation of the t-J model is obtained as 



<hj> 



1 



h.c. 



<i,3> 



Only for "true" spins, the J-term in Eq. (fT^)) represents 
the ordinary magnetic energy. At half-filling, the t-J 
Hamiltonian H [Eq. (TH| ] is reduced to the Heisenberg 
Hamiltonian. 



C. Mean- field approximation 

Within the mean-field approximation, the t-J Hamil- 
tonian H [Eq. (|T4|) ] can be decoupled into a holon part, 
, an d a s P m part, ij¥ F , as 

= ieff^(^i+h.c.)-M$>|/ii, (15) 



<»j> 



H™ F = J cS 22 • Sj - B*g}^S? 



(16) 



<hj> 



iT I T^f+^;) ) J cff = J{(i-<5) 2 - 



with i c ff = i(Si • S j -\ 1 

(/> 2 } + 2i</>, and B z s — — 4i(/>, where (■ ■ • ) is an average 
over the ensemble. The hole doping and the particle- 
hole order parameter are defined as 5 = (h\hi) and <fi = 
(h\hi +n ) with r\ = ±i, ±y, respectively. 

The spin part H^ F [Eq. (fTB)) ] corresponds to the 
Hamiltonian describing the isotropic antiferromagnet in 
a uniform magnetic field with the effective exchange con- 
stant J c ff and the effective magnetic field B z s along z 
axis. We note that the direction of the effective field de- 
pends on the choice of the orientation of the pseudospin. 



nian (fTSf and (fTBf . In Sec. IIII Al we calculate the compo- 
nents of magnetization for the Heisenberg antiferromag- 
net 

H = ^(JSi • S j + D*S*S*) - B^S?, (17) 

<i,j> i 

where J is the exchange constant, D x is the anisotropy 
parameter, and B z is the magnetic field along z axis. 
The Heisenberg model (jTTJ) corresponds to the spin part 
of the mean-field Hamiltonian with 



J - J cff , D X -> 0, S Z -» -B|ff. 



(18) 



Using this correspondence, we calculate the staggered 
magnetization of the mean-field Hamiltonian (|15p and 

dm) m sec. mm 



A. Antiferromagnet in a magnetic field 

The magnetic properties of two-dimensional antiferro- 
magnets have been the focus of many preceding theo- 
retical studie a 15 i 21 ' 22 i 23 ' 24 i 25 ' 26 . In particular, the prop- 
erties of a two-dimensional anisotropic antiferromagnet 
in a transverse field were recently studied with rotat- 
ing frame method 26 and in the non-rotating framed, of 
which difference were discussed by Frobrich and Kunta^. 
In this subsection, we calculate the components of mag- 
netization of the Heisenberg model (fTT|) directly in the 
non-rotating frame as in Ref. [HI. The retarded Green's 
functions are defined as 



G?-(t-t>) = -ie(t-t')([sm,Si(t')}) 

= ({S?;Sr)) ; a = +,-,z , 



(19) 



where 9(t) is the step function. The Fourier transform of 
the Green's functions are denoted by {{Sf ; Sj)) u in the 
energy space. The equations of motion for the Green's 
functions given as Eci. (fT9]) are expressed as 



2(S z )6i ■ 

wG?r(fj)=\ V | +(([S?,ff];ST)) u . (20) 

In order to close the system of equations, we adopt the 
Tyablikov (RPA)-decoupling2£ of higher-order Green's 
functions; 



■ (21) 

Fourier transformations to momentum space are given by 



III. DOPING DEPENDENCE OF STAGGERED 
MAGNETIZATION 

In this section, we consider the doping dependence of 
the staggered magnetization of the mean-field Hamilto- 



(22) 



where subscripts [A, v denote sublattice indices (A or 
B). Each sublattice consists of N/2 lattice sites. Fur- 
thermore, to simplify the calculations, we assume that 
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the magnetic components can be defined as (Sf ) 



"OSS,) 



and (S* ) — (5*4) = mZ ^ ue *° ^ e s y m ~ 



metry of the present case. We now rewrite Eq. (|20| in a 
matrix form, 



uil - 



r o 
o r 



G A 
G b 



A a 
Ab 



(23) 



In Eg. 1(23]). the diagonal 6x6 sub-matrix T is given by 



(24) 
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FIG. 2: Components of magnetization m x , m z as function of 
B z at different temperatures T for D x / J = 0.01 (Solid lines) 
and D x / J = (dotted lines). 



with a = zJto z + B z , b = z(J + D x )m x , c q = 
z(J+ ^£' ;E )7qm z , d q = \zD x ^^m z , and e q = zJ^fqm x , 
where z = 4 is the number of nearest neighbors and 
7 q = 7}(cosq x + cos<7 y ) is the Fourier factor for a square 
lattice: we set the lattice constant to be unity. The sub- 



lattice Green's functions arc defined as 



G„ = [ I : M = A, B , 



- I G 



(25) 



with 



Q+- 

G,„. = | G~ | : /,./'= A. B . i2(i! 

while the vectors A M are given by 

/i = A,B, (27) 



Aa/j^A/x 

Ab m <5b^ 



with 




A M „= | | ; /*,i/ = A,B . (28) 



By solving Eq. (j2"3"]h we obtain the Green's functions with 
the poles 



u = 0, ± wi, ± uj 2 , 

= ^(a + c) 2 + 6 2 - d 2 - e 2 , 

lJ 2 



= v / (a-c) 2 +6 2 -d 2 -e 2 . 



(29) 



Dealing with the pole w = is known to be difficult, in 
that the anticommutator Green's function is required 28 . 
However, we can avoid this difficulty by introducing the 



new expression 

_l 1 /"m x (b + e) — m z (a + c + d) + m z u>i 

Gaa ~ &AA = 

m x (b + e) — m z {a + c + d) — m z iO\ 
1 /m x (b — e) — m z (a — c — d) + m z L02 
m x (b — e) — m z (a — c — d) — m z u)2 



UJ + U2 

(30) 

Since the pole u — is absent in this expression, we 
can use the standard spectral theorem 2 ^ Applying this 
theorem to Eq. (f3T)]) . we derive 



2 



m x (b + e) — m z (a + c + d) f3u)\ 

coth — — 

2wi 2 



m x (b — e) — m z (a — c — d) , /3w 2 

coth 

2w 2 2 



(31) 



In Eq. (|3ip . the sum runs over the first Brillouin zone. 
Moreover, using the fact that the commutator Green's 
function has to be regular at the origin^, i.e., 



lim ujG z ^ = 



(32) 



we obtain 

m x = - 

N- 



(b - e){m x (b + e)-m z (a + c + d)} 
2^| 



(b + e){m x (b — e) — m z (a — c — d)} 



(33) 



From Eqs. f3T]) and |33]) . all components of the mag- 
netization, m x and m z , are determined. As seen in the 
results in Fig. [2] the magnetization m x perpendicular 
to the magnetic field decreases with increasing the field 
strength, and disappears around field strength B z / 'J ' « 4 
atT^Ofor D x /J = 0. 
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B. Doping dependence of staggered magnetization 

In this subsection, we consider the staggered magneti- 
zation of the mean-field Hamiltonian (fTS")) and (fll))) . From 
the transformation ([5]), the staggered magnetization in 
the original representation |(sf}| is transformed into the 
pseudospin representation as 



KOI - K(i - h\hi)sf)\ = \(sf)\ - \(h\hiSf)\ 

= \(Sf)\ 
= M. 



(34) 



Equation ([34| shows that the staggered magnetization 
can be regarded as the component of magnetization per- 
pendicular to the effective magnetic field (to the ori- 
entation of pseudospin) in the pseudospin representa- 
tion. When deriving the relation (|34p. we use the fact 
(h\hiSf) — 0, choosing the orientation of pseudospin 
along z axis. Thus the other components of magneti- 
zation cannot fulfill the similar relation to Eq. (fM)) . For 
example, one can easily show that \{s*}\ -/-> \{Sf)\. 

Instead of starting from the Heisenberg model (|17p. 
we now examine the spin part of the mean-field Hamil- 
tonian ([T5]| . After some algebra similar to the one given 
in Sec. IIII Al we obtain 
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M(b + e) - M(a + c + d) , 0Qi 

coth 

2d>i 2 



M(b -g)- M(q -c-ri) cQth /3w 2 



(35) 



M = 



-y 



(b - e){M(b + e) - M(o + 5 + rf)} 



2£ 2 



(6 + e){M(b - e) - Af(3 - 6 - J)} 



2w| 



+ 



u>i = y(a+ c) 2 +b 2 -<P - e 2 , 
£2 = v/ (S — c) 2 + 6 2 — d 2 — e 2 , 



(36) 

(37) 
(38) 



with M = \(Sf)\, a = zJ eS M + B^, b = zJ eS M, 
c q = zJeS'YqM, dq = 0, and e q = zJ e ff7 q M. In addi- 
tion to Eqs. ([33)) and (f3"6")> . the equations for the holon 
Green's functions are needed to obtain a closed system 
of equations. The holon Green's function are defined as 



gij (t-t') = -i6(t-t')({hi(t),h}(t')}) 
= {{hi-,h})). 



(39) 



From the holon part of the mean- field Hamiltonian (fTS")) . 
we obtain the following expression of the holon Green's 
function; 



f 



(40) 
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FIG. 3: (a) Particle-hole order parameter <j> as function of 
doping 8 at T = (independent of t/J). (b) Staggered mag- 
netization M as function of doping at T = for t/J = 2, 
t/J = 3, and t/J = 5. 



with Sq = ztegfq. Using the spectral theorems^, we also 
find 



27V^ 
q 



1 — tanh 



P{Eq - fl) 



27V^ 7q 
q 



1 — tanh 



2 

/3(e q - M) 



(41) 
(42) 



The set of self-consistency equations (|3"5)) . (|3"6"j) . (|4*T1) . and 
(|32) ar e numerically solved. In Figure [31 we show the 
staggered magnetization and the particle-hole order pa- 
rameter at T — as functions of doping 5. The stag- 
gered magnetization decreases with increasing S, and dis- 
appears at the critical doping S c , which is calculated to 
be 5 C w 0.105 for t/J — 3. In usual mean-field theories, 
since the single occupancy constraint ([2]) is treated only 
on the average, the antiferromagnetic order is generally 
overestimated. However, our results derived without the 
constraint are in good accord with the numerical results 
in Ref. HHllilllM 

In this paper, we choose the pseudospin state as |«t)s- 
If the pseudospin state is chosen as NI)s)/2, the 

direction of the effective magnetic field is aligned along 
the x axis, and the staggered magnetization is mapped 
to the component of magnetization \(Sf )| (perpendicular 
to the effective field). 
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IV. SUMMARY 

In this paper, we have developed the pseudospin rep- 
resentation for the t-J model. By introducing pseu- 
dospins associated with empty sites, we rewrote the t- 
J model in terms of local spin and spinless fermionic 
operators. Using this representation, the t-J Hamilto- 
nian can be described without any constraint. The pseu- 
dospin representation of the t-J Hamiltonian was deter- 
mined by the orientation of the pseudospin state (the 
parameters 9, tp). Within a mean-field level, the spin 
part of the obtained model was regarded as the antifer- 
romagnetic Heisenberg Hamiltonian in a uniform mag- 
netic field with the correspondence B z / J — > B* s /J e s = 
-4t0/(J{(l - 5) 2 - cj) 2 } + 2t(f>). The strength and the 
direction of the effective field are determined by the dop- 
ing S and the orientation of pseudospins, respectively. 



Our method yields reasonable value of the critical dop- 
ing S c , at which the antiferromagnetic long-range order 
disappears, in good agreement with the numerical calcu- 
lationo 16 . 17 i 18 i 19 i 20 . 
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